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Quadrupole Susceptibility and Elastic Softening 
due to a Vacancy in Silicon Crystal 
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We investigate the electronic states around a single vacancy in silicon crystal by using the 
Green's function approach. The triply degenerate vacancy states within the band gap are 
found to be extended over a large distance ~ 20 A from the vacancy site and contribute 
to the reciprocal temperature dependence of the quadrupole susceptibility resulting in the 
elastic softening at low temperture. The Curie constant of the quadrupole susceptibility for the 
trigonal mode {Oyz,Ozx, Oxy) is largely enhanced as compared to that for the tetragonal mode 
{O2, 02)- The obtained results are consistent with the recent ultrasonic experiments in silicon 
crystal down to 20 mK. We also calculate the dipole and octupole susceptibilities and find that 
the octupole susceptibilities are extremely enhannced for a specific mode. 
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1. Introduction 

The recent discovery of low temperature elastic soft- 
ening in silicon (Si) crystal^-* has stimulated much in- 
terest in the investigations not only for physical prop- 
erties but also for industrial applications as the soften- 
ing is closely related to the Si vacancy concentration.^' 
The elastic constants of the both (Cn — Ci2)/2 and C44 
modes show reciprocal temperature dependence below 
20 K down to 20 mK, which implies the existence of 
triply degenerate groundstates which couple to the strain 
caused by the ultrasound.^' The softening of non-doped 
Si is independent of the external magnetic fields up to 10 
T and is attributed to the vacancy with the non-magnetic 
neutral charge state V'^, while that of boron (B)-doped 
Si is suppressed by the external magnetic fields of the 
order of 1 T and is attributed to the vacancy state 
with the valence -1-1 and the spin 1/2.^' 

A huge number of researches on the single vacancy 
in Si crystal have been done in the both experimental'^"^' 
and theoreticaF"^"' viewpoints. Early experiments of the 
electron paramagnetic resonance (EPR)"^' revealed the 
local lattice distortion around the positive charge state 
y+ and negative charge state V~ . Schliiter et al.^' pre- 
dicted that the negative- U effect due to the Jahn- Teller 
distortion results in the 1^++ groundstate with the ex- 
cited state; which is the well-known standard theory 
of the single vacancy in Si crystal and is experimentally 
confirmed by Wat kins and Troxell.^' The negative charge 
state with the Jahn- Teller distortion has also been 
confirmed by the EPR'^' and the electron nuclear dou- 
ble resonance (ENDOR)^' measurements. Recent first- 
principle calculations^^"^''' with supercells up to 1000- 
atoms have shown that the calculated symmetry of the 

*E-mail address; takemi@phys.sc.niigata-u.ac.jp 



relaxed vacancy state depends on computational details 
such as supercell size and k-point sampling. Although the 
detailed symmetry of the charge state is still controversial 
issue, the first-principle calculations have concluded that 
the degeneracies of the charge states are resolved due 
to the Jahn- Teller distortions except the doubly positive 
charge state T/++.1'*' More recently, the elastic proper- 
ties in a Si vacancy have also been discussed by using the 
first-principle calculations.^^' ^^' 

In the EPR experiments,'^' the vacancies were in- 
tentionally created by electron-beam and 7-ray irradia- 
tions, and then, the observed local distortions in the EPR 
are considered to be caused by the cooperative Jahn- 
Teller effect among the highly irradiated vacancies being 
strongly coupled each other. ^' Such vacancies are well 
described by the first-principle calculations with super- 
cells where the vacancy concentration is ^ 1/1000. In 
the ultrasonic experiments, however, the charge state of 
the thermally created vacancies with the extremely low 
concentration in the FZ Si crystal, which is believed less 
than 10^^ cm~'^, shows no sign of the local distortion, 
although a small inter-vacancy coupling ~ 2 K is also 
observed.^' 

When only a single vacancy exists in the infinite Si 
crystal, T^i-point symmetry should be preserved against 
the Jahn- Teller distortions and the 3-fold orbital degen- 
eracy should remain in the vacancy groundstates. Even 
in the case with a finite vacancy concentration, such a 
single vacancy state without the Jahn- Teller distortions 
is expected to be a good starting point to describe the 
properties for finite temperature, where the cooperative 
Jahn- Teller effects due to the small inter-vacancy cou- 
pling are not significant. Such situation is realized in the 
FZ Si crystal above at least 20 mK and are observed 
in the ultrasonic experiments. Therefore we need a the- 
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oretical study for a single vacancy in infinite Si crystal 
without the Jahn- Teller distortions to describe the low 
temperature elastic softening in the FZ Si crystal. 

The similar elastic softenings have been widely ob- 
served in the 4/ electron systems/^' where the Curie- 
like behavior of the quadrupole susceptibilities due to the 
orbital (quadrupole) degeneracy of the 4/ electrons are 
responsible for the elastic softenings via the quadrupole- 
strain coupling. Similarly to the 4/ electron systems, 
the orbital degeneracy of the localized vacancy states 
in the Si crystal is expected to result in the low tem- 
perature softenings observed in the ultrasonic experi- 
ments. Recently, Matsuura and Miyake^^^ have studied 
the quadrupole susceptibility due to a Si vacancy on the 
basis of a cluster model for the dangling-bond orbitals 
in the Si vacancy. "'^^ Yamakawa et al.^^'^^^ have studied 
the effect of the coupling between the electrons of the 
dangling-bonds and the Jahn- Teller distortions by using 
the exact diagonalization for the cluster model developed 
by Schliiter et al.^' and have revealed that the degener- 
acy of the vacancy state remains against the Jahn- Teller 
effect because of the strong quantum fluctuation due to 
the non-adiabatic effect, in contrast to the case with the 
adiabatic approximation.^' However, the effect of the 
spatial extension of the vacancy state, which is important 
to determine the absolute value of the elastic softening, 
has not been discussed in above cluster models. The ex- 
perimental group^-* has claimed that the widely extended 
vacancy state with the effective radius^^^ a ~ 5 A is re- 
sponsible for the huge quadrupole susceptibility in pro- 
portion to which contributes to the elastic constant 
via the quadrupole-strain coupling. 

In this study, we focus on the charge state which 
has non-magnetic groundstate and is realized in the non- 
doped Si as confirmed by the ultrasonic experiment. In 
this case, the spin-orbit interaction, which is crucial for 
the magnetic groundstates in the and states, is 
irrelevant. In the groundstate of V^, the two electrons 
are occupied in the 3-fold degenerate orbital states with 
intra-orbital and inter-orbital spin-singlet configurations 
which are responsible for the electric multipole suscepti- 
bilities such as the quadrupole and octupole. 

The purpose of this paper is to clarify the effect of the 
spatial extension of the vacancy state on the quadrupole 
susceptibility which contributes to the elastic constant 
at low temperature. For this purpose, we determine the 
electronic state around a single vacancy in the infinite 
Si crystal by using the Green's function approach. ^^'^'^^ 
By virtue of this approach, we can discuss the electronic 
state up to 30 A from the vacancy site, where more than 
4000 Si atoms are included. Using the obtained Green's 
functions, we calculate the quadrupole susceptibility on 
the basis of the linear response theory. 



2. Model and Formulation 

The model Hamiltonian consists of the tight-binding 
Hamiltonian Hq and a vacancy potential given by 



H = Hn 
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and c|.j^j is that for 
8). In eq. (2.2), 



where c]^ is a creation operator for an electron at site 

i and orbital a (= s, Px, Py, Pz), "^'^ 
wave vector k and band m (= 1 
the tight-binding parameters t"^^ are written by the 
Slatcr-Koster parameters and determined so as to fit the 
tight-binding band energies ekm to the LDA band ener- 
gies^'*-' for the Si crystal as shown in Fig. 1. The explicit 
values of the parameters are as follows: atomic levels 
eg = —4.778 eV, ep = 1.218 eV, Ist-neighbor hopping 
integrals 4si — —2.104 eV, 4pi = —1.788 eV, ippcn — 
-2.810 eV, ipp^i = -0.743 eV, 2nd-neighbor hopping 
integrals, which gives hopping between the dangling 
bonds, 4s2 = 0.092 eV, 4p2 = 0.112 eV, tpp„2 = 
—0.386 eV, ippT2 = —0.116 eV. The vacancy potential 
i?v excludes electrons from the vacancy site by raising 
the energy levels A for the orbitals belong to the vacancy 
site. For A — > oo, no electron exists at the vacancy site 
and then an effective vacancy state is realized. 
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Fig. 1. (Color online) The band structure for the Si crystal calcu- 
lated from the LDA^'*' (closed circles) and from the tight-binding 
model (solid lines). 



In the absence of the vacancy (A = 0), the Green's 
function for the perfect crystal is described as 

„ u„„(k)M^„(k) 



km 
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where Uam (k) is the eigenvector for the energy band ekm 
and orbital a given in eq. (2.2). In the presence of the 
vacancy (A ^ 0), the Green's function is obtained by 
solving the Dyson's equations which can be written in 
the 4x4 matrix representation as 



(2.5) 

with the vacancy potential matrix (A)^^ = A(5c^, where 
(G° )q/3 = G""^ is the Green's function for A = given 

in eq. (2.4) and {Gij)ap = C^^ is the corresponding 
Green's function for A 7^ 0. In the limit A — > 00, Gy 
with 2 = and/or j — 0, and then 



GESoAGo, ^ 0. (2.6) 
By using eq. (2.6) in eq. (2.5), Gij is obtained by G° as 

1O ^ -1 



G 
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According to the Lchmann representation, G'^^ is de- 
scribed by using the spectral function and the ex- 
citation energy Ei^ as 



A 



El' 



(2.; 



Results 



3. 1 Density of States 

Now, we calculate g'^"^ 
the k summation with 20 x 20 x 10 = 4000 mesh points, 
and substitute it into eq. (2.7) to obtain G"^^ 



y in eq. (2.4) by performing 



which yields 



^^^^ ^Q- (2-8)- By using the obtained val- 

ues of Al^ji^ and Ei, we calculate the local density of 
states (DOS) at site i in the presence of the vacancy, 
P^{^) = E/a ^^Lia^i^-Ei)- Fig. 2 shows p^{uj) together 
with the DOS for the perfect crystal without vacancy, 
where the top of the valence band is set to the energy 
origin. We find two remarkable localized levels : one sits 
in the band gap and the other sits in the valence band. 
By summing the contribution to each level from all sites, 
the total weight of each state is 3 for the former level 
and 1 for the latter level. Therefore, we find that the for- 
mer and the latter levels correspond to T2 triplet states 
with the energy Et2 = 0.44 eV and Ai singlet state with 
the energy Ea^ = —0.12 eV, respectively. These localized 
levels are occupied by 4 electrons in the V° state and by 3 
electrons in the state, respectively. In the both cases, 
the chemical potential n is close to the T2 level at low 
temperature. Then in the following, we focus on the T2 
triplet states which exclusively contribute to thermody- 
namic quantities such as the quadrupole susceptibilities 
at low temperature. 

The inset of Fig. 3 shows the local DOS for the T2 
triplet states at site i, N^^ — J2a ^aion ^ function of 
the distance Ri from the vacancy site. We can see that, 
with increasing i?, Nip^ exponentially decreases with de- 
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Fig. 2. (Color online) The local density of states in the presence 
of the vacancy (solid lines) together with those for the perfect 
crystal without vacancy (broken lines) at 1st (top), 2nd (middle) 
and 3rd (bottom) neighbor Si sites from the vacancy site. 



cay length ^ of several A accompanied by a complicated 
oscillation. We also calculate the integrated local DOS for 
T2 triplet states up to the radius Rc from the vacancy 
site, Nt^ = "Y^f" N'^p^ , which is plotted as a function of 
Rc in Fig. 3. When Rc ^ 00, we can see ^ 3 as ex- 
pected. We note that a large part (~ 80 %) of T2 triplet 
states is concentrated on the region i? < 5 A with 34 
atoms, but a small part (~ 20 7.) of the widely extended 
vacancy states i? = 5 ~ 20 A with ~ 10'^ atoms play 
crucial roles for thermodynamic quantities such as the 
quadrupole susceptibilities as will be shown later. 

3.2 Quadrupole susceptibility 

The low temperature elastic softening is considered 
to be caused by the interaction between the electric 
quadrupole and the elastic strain, which is explicitly 
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Fig. 3. (Color online) The integrated local DOS for T2 triplet 
states up to the radius Rc from the vacancy site. The inset shows 
the local DOS for the T2 triplet states as a function of the dis- 
tance R from the vacancy site. The solid, dot and dashed lines 
are functions exp(— i?/^) with ^ = 2.4, ^ = 2.0 and ^ = 1.6, 
respectively. 



given by 



V 



grTrer{t) 



(3.1) 



with a quadrupole-strain coupling constant gr- Here, 
er(i) is a strain with mode T excited by the ultrasonic 
wave and couples to a quadrupole operator, 



Tr = e^^(m|0r(r)|j/3)c|„Cj-^, 

ij a (3 
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with (r) , given by 
(3z2 
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for 




0r(r) 



where the electric quadrupoles of O2 , 0| couple to the 
tetragonal strain e„,e„ and Oyz,Ozx,Oxy couple to the 
trigonal strain eyz, ezx, £xy, respectively. Since the radius 
of the atomic orbital of each Si atom at site i is smaller 
than the distance Ri from the vacancy site to the Si 
atom, we can approximate (za|0r(r)|j/3) ~ 0r(ri)(Sij(SQ/3, 
and then eq. (3.1) is rewritten as 



HQs{t) ^ -g^ee^{t)y2Mr^)cl 



(3.3) 



The relationship between the elastic constant Cr (T) 
and the quadrupole susceptibility xr{T) at temperature 
T is given in the second order perturbation w.r.t. gr'. 

Cr{T) = C° - N^glxriT), (3.4) 



where Cp is the background of the elastic constant 
and is the number of vacancies. We note that 
the quadrupole susceptibilities of O2 , 0\ contribute to 
the elastic constant of (C\\ — Ci2)/'2,, while those of 
Oyz,Ozx,Oxy contribute to C44. In the linear response 
theory, the quadrupole susceptibility is given by 



Xr(T) = Re / 
"'0 



dt 
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with the fermi distribution function f{E) = 1/(6^*^ ■^^^' + 
1), where Xt{T) consists of the Curie term [Ei = Ei>) and 
the Van-Vleck term [Ei ^ Eii ). As the Van-Vleck term is 
almost T-independent at low temperature, we focus only 
on the Curie term due to the degenerate T2 triplet states 
{El = El' = Ex-i), which is explicitly given by 

E ^r(rO0r(r,)^S.,/-^ /(^Tj/(-i^Tj 



Xr(T) 



j0ia ' 



T 



(3.5) 



with F{n) — n{6 — n)/36, where n — 6/(£'t'2) is the 
occupation number of electrons in the T2 triplet states, 
and n = 2 for the state. Substituting eq. (3.5) into 
eq. (3.4), we obtain the elastic constant which shows the 
reciprocal T dependence, Cr(T) oc —NvgYKrF{n)/T, 
where the absolute value of the elastic softening is deter- 
mined by the Curie constant of the quadrupole suscepti- 
bility given by 



Kr 



TE^r(r. 

ij ap 



(3.6) 



By using eq. (3.6), we calculate with the site 
summation («, j) up to the distance Rc from the vacancy 
site as done for the calculation of , where Kr consists 
of the site diagonal contribution (i = j) and the site off- 
diagonal contribution (i 7^ j). Fig. 4 shows together 
with the site diagonal and site off-diagonal contributions 
for Kr as functions of Rc- To estimate the Rc 00 
extrapolated value of Kp, we assume a fitting function 

KriRc) « i^r(oo)[l-{i??/(2e,') + i?c/er + l}e-^=/«-].^^) 
By comparing the calculated value of Kr with the fit- 
ting function as shown in Fig. 4, we obtain the ex- 
trapolated values: Kr{oo) — 23.28 e^A"' (trigonal) and 
Kr{oo) = 2.51 e^A^ (tetragonal), together with the de- 
cay length: = 3.77 (trigonal) and ~ 3.56 A (tetrag- 
onal). We find that A'r for the trigonal mode is about 10 
times larger than that for the tetragonal mode as shown 
in Fig. 4, where the site diagonal contribution for the 
trigonal mode is almost the same as that for the tetrago- 
nal mode, while the site off-diagonal contribution for the 
trigonal mode is much larger than that for the tetrago- 
nal mode. The result is consistent with the experimental 
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Fig. 4. (Color online) The Curie constant Kr for the trigonal 
mode {Oyz,Ozx,Oxy) (a) and the tetragonal mode (OjjOg) (b) 
as functions of Rc- The open circles, dotted and dashed lines 
correspond to the total, site diagonal and site ofE-diagonal con- 
tributions, respectively. The solid lines show fitting functions for 
the total contributions (see in the text). 



Fig. 5. (Color online) The Curie constants Ky for various mul- 
tipole susceptibilities as functions of Rc- The closed circles (•), 
open circles (o), squares (□), diamonds (o) and triangles (a) 
represent the Kp for the dipole Tx, quadrupole Oxy, octupole 
Txyz,T°' and r*^ susceptibilities, respectively. 



result, where the softening of C44 is considerably larger 
than that of (Cn — Ci2)/2,^) if the coupling constants 
(7r for the both modes, which have not been determined 
so far, are of the same order of magnitude. 

Remarkably, the widely extended vacancy states on 
the region i? = 5 - 20 A contribute to about 90 7. of the 
absolute value of Kr as shown in Fig. 4, although they 
contribute to only 20 7, of the local DOS as shown in Fig. 
3. In fact, the decay length for the Curie constant Kr 
is about twice larger than the decay length ^ for the 
local DOS (see also the inset of Fig. 3). We note that the 
effective quadrupole moments estimated by {KrY^^ are 
4.82 eA^ (trigonal) and 1.58 eA^ (tetragonal), which are 
considerably larger than molecular quadrupole moments 
of the order of 0.1 eA^.^e) 

3.3 Multipole susceptibility 

It is expected that the widely extended vacancy 
states also yield the extreme enhancement of the other 
multipole susceptibilities. Then, we also calculate the 
electric dipole and octupole susceptibilities by using the 
same method as the quadrupole susceptibilities. The 
Curie constants Kr for the corresponding multipole sus- 
ceptibilities are obtained by replacing 0r(r) in eq. (3.2) 
with 

{X for Tx 
y for Ty 

z for 
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for the electric octupole. 

Fig. 5 shows Kr for several multipole susceptibili- 
ties as functions of Rc. It is found that Kr for one or- 
der higher multipole susceptibilities show about 10 times 
larger enhancement. In the octupole susceptibilities, the 
site diagonal contributions for r^yz ,t" are the same 
order of 10 e^A^ (not shown), while the site off-diagonal 
contribution for r" is much larger than that for the other 
mode resulting in the extreme enhancement of r". We 
note that Kr for r'^ decreases with increasing Rc be- 
cause its site off-diagonal contribution is negative. 

4. Summary and Discussion 

In summary, we have investigated the electronic state 
around a single vacancy in infinite Si crystal on the ba- 
sis of the Green's function approach. It has been found 
that the T2 triplet vacancy states within the band gap 
are widely extended up to 20 A and are responsible for 
the extreme enhancement of the Curie constant of the 
quadrupole susceptibilities resulting in the elastic soft- 
ening at low temperature. The Curie constant for the 
trigonal mode {Oyz, Ozxi Oxy) is considerably larger than 
that for the tetragonal mode (O21 0|). These results are 
consistent with the low temperature elastic softening ob- 
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served in the ultrasonic experiments.^' We have also cal- 
culated the other multipole susceptibilities and found 
that the remarkable enhancement of the Curie constant 
for the octupole susceptibilities especially in r" mode; 
which is expected to be observed in future experiments. 

In the present study, wc have assumed that the T^- 
point symmetry remaining the orbital degeneracy is pre- 
served in a single vacancy in infinite Si crystal. Even 
in this lattice relaxation with keeping the sym- 

metry might take place; which results in a modification 
of the tight-binding parameters around the vacancy. To 
discuss this effect, we need a first-principle calculation 
with keeping the symmetry. A recent molecular dynam- 
ics simulation combined with a first-principle calculation 
revealed that such a high symmetric vacancy state is re- 
alized at a finite temperature. ^^-^ The explicit calculation 
of the quadrupole susceptibility including the effect of 
the lattice relaxation is an important future problem. In 
addition, the explicit estimation of the quadrupole-strain 
coupling gr is also important to determine the absolute 
value of the elastic softening. 

At a finite vacancy concentration, the effect of the 
inter- vacancy interaction is considered to be also impor- 
tant. In fact, the ultrasonic experiments revealed that 
there exists the antiferro-type inter-vacancy interaction 
of order of 2 K.^' On the basis of the Green's function 
approach employed in the present study, we can discuss 
the effect of the inter-vacancy interaction by introduc- 
ing the spatially separated two vacancies in infinite Si 
crystal. Such a calculation is under the way. 

The effects of the Coulomb interaction and the cou- 
pling between electrons and Jahn- Teller phonons, which 
have not hccn considered in the present study, are cru- 
cial to deteririine the many-body groundstate in a Si va- 
cancy and have been intensively discussed from the clus- 
ter model calculations. Actually, the non-magnetic 
groundstate in V'^ with 3-fold orbital degeneracy has 
been obtained due to the effect of nonadiabatic couplings 
between electrons and Jahn- Teller phonons,^"' in con- 
trast to the present study where the spin degeneracy is 
unresolved as an artifact of the noninteracting calcula- 
tion. In addition, the effect of the spin-orbit interaction 
neglected in this study is also crucial for the magnetic 
groundstates of and V~ which are expected to be 
realized in the B-doped^^' and P-doped Si, respectively. 
On the basis of the Green's functions for a Si vacancy 
obtained by the present study, we can discuss the many- 
body effect by including the selfenergy corrections due 
to the Coulomb interaction and the electron-phonon cou- 
pling. The detailed results of the selfenergy corrections in 
the presence of the spin-orbit interaction will be reported 
in a subsequence paper. 
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into eq.(3.6) and replacing the i-summation with the rj- 
integral up to |rj| = iic, we can obtain the fitting function 
of Kr{Rc). 
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